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Abstract 



O 

The masses of uudds and uuddd pentaquarks are evaluated in a framework of the 
Q , Effective Hamiltonian approach to QCD using the Jaffe-Wilczek [ud] 2 q approxi- 

mation. The mass of the [uo?] 2 s state is found to be ~ 400 MeV higher than the 
observed O + (1540) mass. 

m 

> 

Recently several collaborations reported the observation of a very narrow peak in the 
K + N invariant mass distribution The reported masses have been clustered around 

1540 MeV, with widths less than 25 MeV. These results are in remarkable agreement with 
^ ! a chiral soliton model prediction [Z| of a G + state at 1530 MeV and a width less than 15 
MeV \ 

i-Q ' 

In the sense of the quark model such a state is manifestly exotic - it is not a simple three- 
quark state but rather the pentaquark uudds state. If so, the discovery of 6 + provides 
an opportunity to refine our understanding of nonperturbative quark dynamics at low 
energy. The B-hyperon has hypercharge Y = 2 and third component of isospin I3 = 0. 
The apparent absence of the ^3 = +1, ++ in K + p argues against 1 = 1, therefore it is 
usually assumed the 9 to be an isosinglet, although other suggestions have been made in 
the literature [TQ| . The other quantum numbers are not established yet. The uncorrelated 
quark models, in which all quarks are in the same spatial wave function, naturally predict 
the ground state energy of a J p = \ pentaquark to be lower than that of a J p = | + 
one. Several suggestions were made to reverse this order jTTj, [T2] but no quantitative 
microscopic evaluations have been performed yet. The QCD sum rules predict a negative 
parity G + of mass ~ 1.5 GeV, while no positive parity state was found [T^j. The lattice 
QCD study also predicts that the parity of the 9(1540) is most likely negative [Tlj.fTT 



Jaffe and Wilczek proposed jTTj that for 0(1540) with I p = + and other q A q baryons 
the four quarks are bound into two scalar, singlet isospin diquarks combined with the 
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1 The problem of the 8 + width is crucial since the widths beyond the few-MeV level seem to be 
excluded by the K + N phase-shift analysis [B| and total cross sections jHj- 
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Table 1: Comparison of the EH approach predictions for spin-averaged masses of nucle- 
ons, hyperons, and pentaquark with lattice results and experimental data (assuming the 
@ + quantum numbers are I p = + ). The masses are given in units of GeV. 
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antiquark into a color singlet. The problem of diquarks has a long history and sound 
physical motivation ^B]. Both gluon and pion exchanges in the singlet spin and isospin 
ud state result in the strong short-range attraction which may keep u and d quarks close 
together thus forming almost point-like diquark. 

The purpose of this letter is to test this interpretation quantitatively using the Effective 
Hamiltonian (EH) approach in QCD An attractive feature of this formalism is that 
it contains the minimal number of input parameters: current (or pole) quark masses, 
the string tension a and the strong coupling constant a s , and does not contain fitting 
parameters as e.g. the total subtraction constant in the Hamiltonian. Therefore the 
EH method is the best way to compute the masses of yet unknown quark systems like 
pentaquarks. 

In order to illustrate the accuracy of the method, we quote in Table [T] a few EH results 
for baryons compared with the lattice ones and experiment. One can observe that the 
accuracy of the EH method for the three-quark systems is ~ 100 MeV or better. We 
expect the same accuracy for the diquark-diquark-(anti) quark system considered in [TT] . 

The EH for the three constituents has the form 



where Hq is the kinetic energy operator, V is the sum of the perturbative one-gluon 
exchange potentials and the string potential which is proportional to the total length 
of the string, i.e to the sum of the distances of (anti)quark or diquarks from the string 
junction point. The dynamical masses fa are expressed in terms of the current masses 
m i from the condition of the minimum of the hadron mass M$ as function of fa 3 : 

dMji\mi, fa, n u mi \ • / n \ L . ,,, , ^) 



dfa 

E (fa) being eigenvalue of the operator H Q + V. The physical mass M H of a hadron is 

M H = MP + Ys C " ( 3 ) 



2 This work 

3 Technically, this is done using the auxiliary field approach to get rid of the square root term in 
the Lagrangian [23, |22| . Applied to the QCD Lagrangian, this technique yields the EH for mesons 
or baryons depending on auxiliary fields fa. In practice, these fields are finally treated as c- numbers 
determined from (J2J. 



2 



a 



b 



Figure 1: The Jaffe-Wilczek reduction of the uuddq pentaquark (a) to the 
effective [ud] 2 q problem (b) 



where the constants Q have the meaning of the constituents self energies and are ex- 
pressed in terms of string tension a 23J. The self-energy is due to constituent spin 
interaction with the vacuum background fields and equals zero for any scalar diquark. 

In the quark model five quarks are connected by seven strings (Fig. la). In the diquark 
approximation the short legs on this figure shrink to points and the five-quark system 
effectively reduces to the three-body one (Fig. lb), studied within the EH approach 
in [21], Consider a pentaquark consisting of two identical diquarks with current 
mass rri[ u d] and antiquark with current mass mq (q = d,s). The three-body problem 
is conveniently solved in the hyperspherical formalism j2E]- The wave function ip(p, A) 
expressed in terms of the Jacobi coordinates p and A can be written in a symbolical 
shorthand as 

iP(p,X) = Y^^(R)Y [K] (n). (4) 

K 

where Y\}o are eigen functions (the hyperspherical harmonics) [2Zj of the angular momen- 
tum operator K(Q) on the 6-dimensional sphere: K 2 (Q)Y[k] = —KiK + 4)Yr#n, with K 
being the grand orbital momentum. For identical diquarks, like [ud] 2 , the lightest state 
must have a wave function antisymmetric under diquark space exchange. There are two 
possible pentaquark wave functions antisymmetric under diquark exchange, one (with 
lower energy) corresponding to the total orbital momentum L = 1, and the second one 
(with higher energy) corresponding to L = 0. For a state with L = 1, l p = 1. l\ = the 
wave function in the lowest hyperspherical approximation K = 1 reads 



i> = R- 5 / 2 Xl (R) Ul (n), Ul (Q) = J-sm6-(n p ) z 

V 7T 



(5) 



where R 2 = p 2 + A 2 . For a state with L = 0, l p — 1, l\ = 1 the wave function in the 
lowest hyperspherical approximation K = 2 is 

i/j = R~ 5 / 2 X 2(R)u 2 (n), u 2 (Q) = ^ sin 6 cos 9 ■ (n p n x ). (6) 



The Schrodinger equation written in terms of the variable x = ^/JiR, where /i is an 
arbitrary scale of mass dimension which drops off in the final expressions, reads: 



d 2 XK{x) , n \„ , a K L (K+I)(K , , 



dx 2 



+ 2 



E Q H 



x ' 2a; 
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Xi^(x) = 0., (7) 



with the boundary condition Xk(%) ~ (9(s 5 / 2+K ) as a; — > and the asymptotic behavior 
XJf(s) ~ Ai((2b K ) 1/3 x) as x -> oo. In Eq. (0) 



a K = RtJJi ■ J Vc ou i omb (ri, r 2 , r 3 ) • • dO, 
&x = • y Vrtring(ri, r 2 , r 3 ) • it^ • dfi, 



where 

2 1 

Vcouiomb(T , i, r a , r 3 ) = --a, • — , (9) 

K] J 

and explicit expression of V stT i ng (ri, r 2 , r 3 ) in terms of Jacobi variables is given in |T8*j . 
Note, that in the Y-shape in Fig. 16, the strings meet at 120° in order to insure the 
minimum energy. This shape moves continuously to a two-legs configuration where the 
legs meet at an angle larger than 120°. 

The mass of the + obviously depends on m\ U( n and m s . The current masses of the 
light quarks are relatively well-known: ui Ui & ~ 0, m s w 170 MeV. The effective mass 
of the diquark rri[ u d} is basically unknown. In principle, this mass could be computed 
dynamically. Instead, one can estimate from the nucleon mass calculations in 

the quark-diquark approximation. In this case, neglecting for simplicity Coulomb-like 
interaction, one obtains for the mass of the nucleon: 

M N = MP - ^, (10) 

where M®' and \i& are defined from the minimum condition, 

um N _ um N _ 



In Eq. ffTTj) 



mP = ^ + ^±^ + E 0i E = 1 (^.^±^) 1/3 (12) 
2/i[«d] 2 \2 H{ud]tJ-d J 

with 7 = 2.338 being the first zero of the Airy function: Ai(— 7) = 0. Equating Mn to 
the experimental value of the N — A center of gravity (M^ = 1.085 GeV), we obtain 
that m\ u d] varies in the interval < mry < 300 MeV, when a varies in the interval 
0.15 GeV 2 < a < 0.17 GeV 2 . The last value of a is preferred by meson spectroscopy, 
while the former one follows from the lattice calculations of Ref. j2B|- The lowest 
uudds pentaquark corresponds to the case when one fixes the effective mass of a di- 
quark m[ u d} = 0. In what follows, we use = 0, a = 0.15 GeV 2 , but explicitly include 
the Coulomb-like interaction between quark and diquarks. 



4 



Table 2: The pentaquark results in the quark-diquark-diquark approximations. Shown 
are the constituent masses //[ uc q, fi q , eigenvalues E in Eq. (J2J), and masses of [udfq states 
(q = d, s) for J p = | . = 0, m d = 0, m s = 0.17 GeV 
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[udfd \~ 
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The calculated masses of [ttrf] 2 s pentaquarks are given in Table El For illustration we 
show also the masses of [udfd pentaquarks which have been identified in pi] with the 
otherwise perplexing Roper resonance. The states with K = 2 are always higher than 
those with K = 1 by ~ 300 MeV, the difference being mainly due to the kinetic energy 
term in (J7J). Increasing a s up to 0.6 (the value used in the Capstick-Isgur model [2H]) 
decreases the [mc/] 2 s mass by ~ 100 MeV. The hyperfine interaction due to rj exchange 
between diquarks and strange antiquark can lower the G + energy by ~ 100 MeV or less 
[12]. Recall that the accuracy of our approach for pentaquarks is expected to be ~ 100 
MeV. Since the calculated value is the lower bound on the pentaquark mass, we conclude 
that the lightest uudds state is still ~ 400 MeV higher than the observed G + (1540) state. 
Note that [udfd pentaquarks lie ~ 100 MeV lower than [«(i] 2 s pentaquarks. This is 
the consequence of the different current masses m rf - and m s . Also note, that for the \ 
pentaquark the [ltd] [tisj-u configuration with K min = is more preferable giving the mass 
> 1.8 GeV. 

In conclusions, we have presented the first dynamical calculations of the pentaquark 
masses within the Jaffe-Wilczek approximation. Our predictions are universal in the 
sense that they use only two parameters inherent to QCD: the string tension a and the 
strong coupling constant a s . These results do not require any additive mass shift used 
in constituent quark models. Our predictions can be viewed as lower bounds for the 
masses of pentaquark states. We showed that the + mass is much higher than the 
observed one. This implies that O + (1540) can not be explained in terms of the standard 
string interaction between quarks in the diquark-diquark-antiquark picture. Note that 
the chiral interactions (e.g. pion exchanges) are automatically taken into account inside 
diquarks, while are not possible between diquarks and strange antiquark. Therefore 
a drastic modification of the present results requires a completely new dynamics, ei- 
ther the chiral soliton type, as in [2], or other quark clusters, like q A — q [30] or q 2 s — q 2 |31j . 

The authors are grateful to K.G.Boreskov for careful reading of the manuscript and 
valuable suggestions. This work was supported by RFBR grant # 03-02-17345 and PRF 
grant for leading scientific schools # 1774.2003.2. 
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